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ABSTRACT 

We analyze the gamma-ray halo around stacked AGNs reported in Ap.J.Lett., 2010, 722, L39. First, we show that the angular dis- 
tribution of y-rays around the stacked AGNs is consistent with the angular distribution of the y-rays around the Crab pulsar, which 
is a point source for Fermi/LAT. This makes it unlikely that the halo is caused by an electromagnetic cascade of TeV photons in the 
intergalactic space. We then compare the angular distribution of y-rays around the stacked AGNs with the point-spread function (PSF) 
of FermifLKT and confirm the existence of an excess above the PSF. However, we demonstrate that the magnitude and the angular size 
of this elfect is different for photons converted in the front and back parts of the Fermi/LAT instrument, and thus is an instrumental 
effect. 
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High-energy gamma rays propagating through the 
Universe produce electron-positron p airs in interactions 
with the extragalactic b a ckground light ( Kneiske et al.L l2004t 
iFranceschini et al.L 120081: IStecker & ScuUvi |2009|) . The leading 
particle of the pair then upscatters the background photons 
through the inverse Compton effect, creating an electromagnetic 
cascade. If there is a weak intergalactic magnetic field, electrons 
and positrons in the cascade are deflected from the original di- 
rection before being converted back to photons. These secondary 
photons create a "halo" of softer gamm a-ray photons around an 
extr agalactic TeV gamma-ray so urce (Aharonianet al., 1994'; 
PIaga..l995:iNeronov & Semikoa ,2007; Takahashi et aU200 



Elviv et aI.L (20091: iDoIag et aI.LT2009l: iNeronov et aI.L 1201 
This process may be used to constrain the parameters of the 
extragalactic magnetic fields ( Plaga, 1995; Neronov & S emikozl 
I2007t iTakahashi et al.L I200I lEIviv et all 120091) . The absorption 
of TeV y-rays leads to the cascade emission in the GeV energy 
range, which is accessible for obse rvations with the Fermi tele- 
scope eronov & Semikozl l2009l) . Non-detection of extended 
emission around blazars with hard intrinsic spectra extending 
to the muIti-TeV energy range by Fermi was used to derive a 
lower boundary on magnetic fields in the intergalactic medium 
at the level of 10"'^ G to 10"'^ G, depending on the assumptions 
about the intr insic spectral properties of the analyzed sources 
( Neronov & VovL .2010: Tavecchio et al. ,2010a.b.: Dolag et al.L 
I2010h . 



lAndo & Kusenkd (1201 Oh recently claimed to posess evi- 
dence for the gamma-ray halos around AGNs from the Fermi 
source catalog. The halo was found in the stacked signal ofAGNs 
selected by certain criteria. It was interpreted to be caused by 
cascading and deflections of TeV photons in the extragalactic 
magnetic fields. Unfortunately, we find that the proposed inter- 
pretation of the halo is unlikely, but a more probable cause is a 



tail in the point-spread function (PSF) of photons pair converted 
in the back thick layer of the Fermi detector 
Our argument consists of two parts. 

First, we compare the stacked AGN signal to the signal of the 
Crab pulsar. This is a bright galactic gamma-ray source whose 
signal in the Fermi energy band consists of two contributions: 
emission from the pulsar and from the associated pulsar wind 
nebula (TWN). The Crab PWN has an angular size ^ 0.05° 
(lHesteiLl2008l) . which is below the angular resolution of the LAT 
telescope onboard Fermi. This means that the Crab PWN is a 
point source for LAT. In addition, this PWN is situated in the 
Galactic anti-center region, where the density of sources in the 
muIti-GeV energy band is low and the diffuse Galactic y-ray 
background is relatively low as well. This means that Crab sig- 
nal is not contaminated by the nearby point sources or by strong 
inhomogeneities of the diffuse background, and the source can 
be considered as an isolated point source. 

The LAT PSF depends on the photon energy. Thus, the shape 
of the angular distribution of y-rays around a source depends 
on the source spectrum. The total Crab pulsar -hPWN spectrum 
in the 10-100 GeV band is well describe d by a power law 
dN/dE ~ E-^ with the photon index F ^ 2 dAbdo et aI.Ll2010l) . 
which is the same as the cumulati ve spectrum of the AGNs de- 
rived by lAndo & Kusenkd (l2010l) . This implies that the point- 
source contribution to the AGN signal should have the same an- 
gular shape as the Crab signal. Thus, a difference between the 
angular distributions of photons around AGN and the Crab PWN 
would indicate a halo around AGNs in addition to the point- 
source contribution. 

Fig. [U shows the angular distribution of photons around the 
stacked AGNs (red) and the Crab source (black). In both cases 
the background is subtracted. The shapes of the two signals 
coincide, which means that the entire stacked AGN signal is 
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Fig. 1. Comparison of the background-subtracted angular distributions of 3 - 10 GeV (left) 10 - 100 GeV (right) y-rays around the 
position of the Crab pulsar (black) and around the stacked AGNs (red). Normalization of the distribution of y-rays around AGN is 
scaled to match the y-ray distribution around the Crab pulsar. In the 3-10 GeV band AGN with redshift z < 0.5 are shown in blue. 



well described by a point-source signal, with no additional halo 
contrib ution. In an update of their analysis, lAndo & Kusenkol 
(120101) have claimed that the halo signal appears only in a sub- 
set of AGN at low redshift, z < 0.5 and only in the 3-10 GeV 
energy band. To allow an explicit comparison with Fig. 4 of 
lAndo & Kusenkol (1201 Ol) . we show in blue the photon distribu- 
tion around AGN at z < 0.5 in the left panel of Fig. [1] which 
corresponds to photon distributions in the 3-10 GeV band. No 
discrepancy between the Crab PWN and AGN profiles is seen in 
this case either. 

Second, we have investigated the nature of the excess in the 
angular distribution of photons around AGN above the LAT PSF 
in the energy range 1 0- 1 00 GeV. y-rays detected by the LAT tele- 
scope are split into two types; photons that are pair-converted 
in the thin front layer of the LAT detector ("front" photons), 
emd photons converted in the thick back layer ("back" photons) 
dAbdo et al.L |2009|) . The shapes of the PSF for these two types 
of photons are significantly different. Taking this into account, 
we spli t the entire photon signal from the set of AGNs consid- 
ered bv lAndo & Kusenkol (12010) into two parts corresponding to 
the front and back converted photons and analyze each part sep- 
arately. We then compare the angular distributions of the front 
and back photons with the corresponding PSFs calculated using 
the DIFFUSE_P6_v3 calibration files by averaging the PSF for 
a given photon energy and incidence angle over the entire set of 
energies and incidence angles of detected photons from the AGN 
set. 

The angular distribution of all (front + back) photons around 
AGN is shown in Fig. |2] together with the corresponding PSF. 
We find an excess in the data at 0.2° <e < 0.9° for the 10-100 
GeV band and in the 0.5° - 1° range in 3-lOGeV band. There is 
a significant difference between our estimate of the background 
level shown in Fig . |2l an d the background level in the Fig. 2 of 



lAndo & Kusenkol ( 1201 Ol) . which is a factor of ^ 3 higher than 
our estimate in the 10-100 GeV band. The constant background 



level shown in Fig. [2] was derived from fitting of the data in the 
range < < 2° around AGN with a model consisting of PSF 
plus a constant. 

We have checked that our estimate of the background level 
is self-consistent in the following way. The total number of 
photons detected by Fermi in the energy band 10-100 GeV in 
the Galactic latitude range \b\ > 10° is A^totai - 5.6 x 10"*. 
Only a small fraction of these photons, A^source - 0.6 x 10"* 
could be associated to the known Fermi sources. The remain- 
ing photons contribute to the diffuse Galactic and extragalac- 
tic backgrounds. Dividing the number of diffuse background 
photons by the solid angle Q spanned by the considered part 
of the sky, one finds the surface brightness of the background 
dN/d6^ - (Motai-A^source)/f^ - 1-5 dcg"^. Multiplying this num- 
ber by the number of AGN considered in the analysis one arrives 
at the background estimate consistent with the one shown in Fig. 
|2] We stress that t his background est i mate i s significantly lower 
than that found bv lAndo & Kusenkol ( 1201 Ol) 

The angular distributions of front and back photons sepa- 
rately are shown in Fig.[3j each with the corresponding PSF. One 
can see that the front photons have more compact PSF, which is 
reasonably compatible with the angular distribution of the front 
photon part of the AGN signal. The fit of the data in the range 
< 6» < 2° by the PSF + background gives x-/dof = 125/98 
{p ~ 3%) in the 10-100 GeV band. Top panels of Fig. 3 show 
the angular distributions of front and back photons in the 3-10 
GeV band. Obviously there is no hint of excess above the PSF 
of the front photons in this energy band. The reduced x~ of the 
fit is ;i'-/dof = 67/72. 

The back photons have a wider PSF, which does not repro- 
duce the back part of the AGN signal in either 3-10 GeV or 10- 
100 GeV energy range. 

Adding front and back photon signals one can verify that 
most of the excess above the DIFFUSE_P6_v3 PSF in Fig |2] is 
caused by the back photons. There is no strong excess above the 
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Fig. 2. Angular distribution of photons around stacked AGNs. Curves show the PSF (thin solid line), the background level (dashed 
line) and PSF plus background (thick solid line). 

PSF in the front photon signal. This i ndicates that the excess o f 
the photons above the PSF, reported in lAndo & Kusenkol ( 1201 Oh . 
cannot be attributed to the extended emission around AGN, be- 
cause in this case an equally strong excess should be detected in 
both the front and back photon signals. 

To quantify the excess in the front and b ack photon signals 
above the DIFFUSE_P6_v3 PSF we followed IXndo & Kusenkol 
(120101) and introduced an additional ^^-Gaussian component 
dN/d9^ ~ exp(-[^?/0haio]'*) in the model of the angular photon 
distribution. This additional component depends on two param- 
eters, /halo (fraction of the source signal contained in the addi- 
tional component) and flhaio (angular size of the core). Fitting 
the angular photon distributions with the improved model we 
found the best-fit values and confidence regions for /haio, ^aio 
separately for the front and back photons. The result is shown in 
Fig. |4] The parameters of the additional components found for 
the front and back photons are clearly different. An addition of 
the new component to the model does not improve the fit of the 
front photon angular distribution in the 3-10 GeV. Only an upper 
bound on the flux in the "halo" component could be derived, as 
is shown in the left panel of Fig.H) This confirms our conclusion 
that the excess above the PSF is not caused by the real extended 
emission around AGN. 

Indeed, the size and the flux of a real halo around AGN 
cannot depend on the way the photons were converted into 
e^e" pairs in the Ferm;/LAT detector. Instead, the excess above 
the PSF in the back photons should be attributed to the im- 
perfect modeling of the real PSF for back photons in the 
DIFFUSE_P6_v3 version of the LAT instrument characteristics. 

To summarize, the angular distribution of photons around the 
stacked AGNs matches well that around the Crab source in both 
the 3-10 GeV and 10 - 100 GeV bands. A separate analysis 
of the front and back converted photons shows that there is no 
evidence for a physical halo around AGN, neither in 3-10 GeV 
nor in 10-100 GeV band. No halo caused by the extragalactic 
magnetic fields is observed. 
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Fig. 3. The same as Fig. 2, but separately for the front (left panel) and back (right panel) photons in 3-10 GeV band (top) and 10-100 
GeV band (bottom). 
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Fig. 4. 68 % (solid) and 99% (dashed) confidence contours for the halo parameters 0haio,/haio derived from the separate analysis of 
front and back photons. 



